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A 


Status from Applications Point of View 


• Applications to Real-World Problems 

- N-S solution of full conf iguration was a big goal in the 80s 

- Numerical procedures and computing hardware have been 
advanced enabling simulation of complex conf igurations 

• Some Examples of Successful Applications 

- Components of liquid rocket engine 

- Hydrodynamics (Submarines, propellers, ...) 

- Cround vehicles (automobile aerodynamics, internal flows...) 

- Biofluid problems (artificial heart, lung, ..,) 

- Some Earth Science problems 

• Current Challenges 

- For integrated systems analysis, computing requirement is very large 
=> Analysis part is still limited to low fidelity approach 

- For high-fidelity analysis, especially involving unsteady flow, long turn-around 
time is often a bottle neck => Acceleration of solution time is required 


Major Drivers of Current Work 
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i To provide computational tools as an economical option for developing 
future space transportation systems (i.e. RLV subsystems development) 

Impact on component design => Rapid turn-around of high-fidelity analysis 
Increase durability /safety ==> Accurate quantification of flow 

(i.e. prediction of low-induced vibration) 

Impact on system performance => More complete systems analysis 

using high-fidelity tools 

> Target 

Turbo-pump component analysis => Entire sub-systems simulation 
Computing requirement is large: 

=>The goal is to achieve 1000 times speed up over what was possible in 1992 
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Objectives 
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• To enhance incompressible flow simulation capability for developing 
aerospace vehicle components, especially, unsteady flow phenomena 
associated with high speed turbo pump. 


advanced 
engine 

phase II + redesign 


engines 


artificial heart 


INS 30 - Incompressible N-S Solver 




** Parallel version : Based on INS3D-UP 

•MPI and MLP parallel versions 

• Structured, overset grid orientation 

• Moving grid capability 

• Based on method of artificial compressibility 

• Both steady-state and time-accurate f ormulations 

• 3 rd and 5 th -order flux difference splitting for convective terms 

• Central differencing for viscous terms 

• One- and two-equations turbulence models 

• Several linear solvers : GMRES, line- relaxation, LU-SGS, 

point relaxation, ILU(O)..,.. 

•HISTORY 

** 1982-1987 Original version of INS3D - Kwak, Chang 
** 1988-1999 Three different versions were devoped : 
INS3D-UP / Rogers, Kiris, Kwak 
INS3D-LU / Yoon, Kwak 
INS3D-FS / Rosenfeld, Kiris, Kwak 
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Time Accurate Formulation 


Time-integration scheme 

Artificial Compressibility Formulation 

* Introduce a pseudo-time level and artificial compressibility 

* Iterate the equations in pseudo-time for each time step until 
incompressibility condition is satisf ied. 


Pressure Projection Method 

* Solve auxiliary velocity field first, then enforce 
incompressibility condition by solving a Poisson equation 
for pressure. 


Artificial Compressibility Method 
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Time-Accurate Formulation 

Discretize the time term in momentum equations using second-order 
three-point backward-difference formula 


dU + dV + dW 
0 £ 07 / dC 


+ dV + dW =Q 3 g~‘-4q-+<r = _„, 


Introduce a pseudo-time level and artificial compressibility. 
Iterate the equations in pseudo-time for each time step until 
incompressibility condition is satisfied. 

^ ^pn + \,m + 1 — — I 
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Time History of Stagnation Point 


Current Challenges 


• Challenges where improvements are needed 

- Time-integration scheme, convergence 

- Moving grid system, zonal connectivity 

- Parallel coding and scalability 

• As the computing resources changed to parallel and distributed 
platforms, computer science aspects become important such as 

- Scalability (algorithmic A implementation) 

- Portability, transparent coding etc. 

• Computing resources 

- "Grid" computing will provide new computing resources for 
problem solving environment 

- High-fidelity flow analysis is likely to be performed using “super 
node" which is largely based on parallel architecture 


Page 5 







Page 6 















Page 8 













Page 9 





Page 10 







Page 1 1 






Page 1 2 







Page 13 







PRESSURE 


VELOCITY MAGNITUDE 


NASA 


- 34.3 Million Points 

- 800 physical time steps in one 
rotation. 

*One physical time-step requires 
less then 20 minutes wall time 
with 80 CPU's on Origin 2000, 

One complete rotation requires 
one-week wall time with 80 CPUs. 
^Currently I/O is through one 
processor. Timing will be improved 
with parallel I/O since time- 
accurate computations are I/O 
intensive. With futher 
improvements several impeller 
rotations can be completed in one 
week. 


Page 14 





Summary 



• Unsteady SSME-rigl start-up procedure from the pump at rest has been 
initiated by using 34.3 Million grid points. 

• Computational model for the SSME-rigl is completed. Moving boundary 
capability is obtained by using DCF module in OVERFLOW-D. 

• MPI /Open MP hybrid parallel code has been benchmarked. 

• MLP shared memory parallelism has been implemented in INS3D, and 
benchmarked. 

• MLP/OpenMP version requires 19-257© less computer time than 
MPI/OpenMP version. Pin-to-node for MLP version is implemented. 407© less 
computer time is required in the new version. 

• Time-accurate features of methods designed for 3-D applications are 
evaluated. An efficient solution procedure is obtained. 

• Work currently underway 

• Unsteady SSME-rigl simulations by using 34.3 Million grid points. 

• Experimental measurements at NASA-MSFC. 
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